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Abstract. A fast method is presented for computing the topological entropy of braids on the 
torus. This work is motivated by the need to analyze large braids when studying two-dimensional 
flows via the braiding of a large number of particle trajectories. Our approach is a generalization 
of Moussafir's technique for braids on the sphere. Previous methods for computing topological 
entropies include the Bestvina-Handel train-track algorithm and matrix representations of the braid 
group. However, the Bestvina-Handel algorithm quickly becomes computationally intractable for 
large braid words, and matrix methods give only lower bounds, which are often poor for large braids. 
Our method is computationally fast and appears to give exponential convergence towards the exact 
entropy. As an illustration we apply our approach to the braiding of both periodic and aperiodic 
trajectories in the sine flow. The efficiency of the method allows us to explore how much extra 
information about flow entropy is encoded in the braid as the number of trajectories becomes large. 
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1. Introduction. Investigation of two-dimensional fluid mixing by topological 
techniques is rapidly gaining popularity [1]. Topological perspectives on mixing either 
involve studying braiding motion of the stirring apparatus itself [5] , or the diagnosis 
of mixing by analyzing braiding of orbits of the flow [10-12, 18]. The quantity that is 
usually of interest is the topological entropy of the braid [4] , which serves as a lower 
bound for the topological entropy of the flow. The topological entropy of the flow is 
related to the exponential growth rate of material lines [17], which has long been a 
favorite measure of mixing quality, though it is by no means the only one [9] . In many 
cases the braid entropy is quite a sharp bound on the flow entropy [8, 10, 12], which 
is one reason why analyzing braids is useful. Another reason is that experimental 
particle trajectory data can be found easily by particle image velocimetry, but it is 
usually very difficult to measure entropies directly from material stretching or by 
computing Lyapunov exponents. 

There are many techniques for calculating braid topological entropies, or lower bounds 
on them, including train-tracks [2,13] , the Burau representation of the braid group [14] , 
and others [15,16]. Where braiding of periodic (or aperiodic) orbits is used to analyze 
a flow one needs to interpret braids that have a both a large number of strands and 
a large number of generators. In this scenario exact methods based on train-tracks 
quickly become prohibitively expensive computationally, and methods based on the 
Burau matrix representation of the braid group usually give very poor lower bounds. 
Recently, however, Moussafir described a fast method for calculating the entropy of a 
braid to an arbitrary precision [16]. This method is based on a Dynnikov coordinate 
representation of a lamination [7]. In this paper we show how to extend Moussafir's 
technique to braids on the torus. We are motivated by the fact that many interest- 
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Fig. 1.1. The spatially periodic domain with n punctures, (a) Any motion of the punctures 
(up to homotopy) can be written as a braid word consisting of a string of o~i , Pi and Ti motions and 
their inverses, (b) It is convenient to define the additional operation cr n , the clockwise exchange of 
the last and first punctures. 



ing dynamical systems are defined in periodic (cylindrical or annular) or biperiodic 
(toroidal) spatial domains. For instance, the alternating sine flow [10] can be analyzed 
from a topological perspective using this approach. 

Figure 11.11 shows the setting for the dynamical system under study. We consider a 
flow on a torus, so that the domain is periodic in both directions. We identify n 
distinguished points that we call punctures. We then consider motions that move 
the punctures in such a way that they always return to their initial configuration, 
possibly having been permuted amongst themselves. The three types of motions that 
we consider are illustrated in Figure fTTTT a): they are 

1. <7i, the clockwise interchange of the ith and (i + l)th puncture; 

2. pi, the ith puncture making a full tour around the vertical periodic direction; 

3. Ti, the ith puncture making a full tour around the horizontal periodic direc- 
tion. 

The inverse of any of these motions is obtained by reversing its direction. The el- 
ementary motions {cri,pi,Ti} are generators of the braid group on n strands on the 
torus [3]. Recognizing the periodicity of the domain, we also define an additional op- 
eration er„ to be the clockwise interchange of the nth puncture with the first puncture. 
To be precise, we mean here the first puncture in the 'copy' of the domain above and 
to the right of the nth puncture, as pictured in Figure rTTT b). so that both periodic 
boundaries are crossed in performing a n . Defining a n in this way keeps both periodic 
directions on an equal footing, and is also convenient in what follows as it is related 
to a translational symmetry for the punctures. 

A sequence of generators, such as Pg 1 <J2P&&&, is called a braid word, and we use the 
convention that the elementary motions in a braid word occur from left to right, so 
that p^ 1 occurs first in our example. By planar braid we mean a braid word that can 
be written using only generators from the set {cti, . . . , cr n -i}, and their inverses. A 
planar braid is equivalent to a braid on the plane with n punctures. In other words, 
a planar braid does not take advantage of the periodic directions. By cylinder braid 
(or annular braid) we mean a braid word that can be written using only generators 
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Fig. 2.1. (a) A representative curve of a lamination on a toroidal domain with five punctures. 
The top and bottom edges, and the left and right edges, are identified so that the picture shows 
exactly ten triangles, (b) The number of crossings with each edge of the triangulation, from which 
the lamination can be reconstructed. 

from either the set {ci, . . . ,a n —x,pi, . . . , p n } or {<7i, . . . , er n _i,Ti, . . . , t„}, and their 
inverses. In other words, a cylinder braid takes advantage of one periodic direction, 
but not the other. Finally, a torus braid is a braid word that is neither a planar braid 
nor a cylinder braid. In this paper we will derive a general method for torus braids, 
which includes planar and cylinder braids as special cases. 

The description of our method here is intended to be accessible to dynamicists and 
requires no specialist understanding of braid groups. The paper is divided as follows. 
In the following section we describe how laminations (equivalence classes of simple 
closed curves) can be encoded by triangulation of the flow domain. Section 03 the 
heart of the paper, gives the details of how this encoding evolves under fundamental 
braiding actions or motions. Several examples are given in Section 0] to illustrate and 
verify the method. In Section [3] we show how to compute the topological entropy 
from evolution of laminations, and demonstrate the rapid convergence. We also look 
at braiding in the sine flow as an example application. We summarize our work in 
Section |H1 and discuss some features and possibilities for further study. 

2. Encoding of laminations by triangulation. We wish to calculate a lower 
bound on how rapidly material lines are stretched in a continuous-time flow based 
on the motion of a finite number of punctures, or a finite set of periodic orbits. 
Our approach, inspired by the method of Moussafir [16], is to study the stretching 
and folding of laminations. A lamination is an equivalence class (under homotopy) of 
simple closed curves that are not homotopic to any part of the boundary (treating the 
punctures as boundaries), and cannot be continuously shrunk to a point. For example, 
a loop that encloses at least two punctures belongs to a lamination. We usually 
represent a lamination by drawing one loop in the equivalence class. An example of 
a lamination on our doubly-periodic toroidal domain is shown in Figure !2"1T a) . 

In order to calculate how a given lamination is transformed under the action of the 
braid, it is necessary to have a way of encoding the lamination, and a method for 
evolving this encoding under the action of the braid. An elegant way of encoding 
the lamination is by triangulating the entire domain and counting the number of 
crossings the lamination makes with the edges of the triangulation. The number of 
crossings for our example lamination is shown in Figure I^TTT b). Note that there are 



4 



M. D. FINN AND J.-L. THIFFEAULT 




Fig. 2.2. Notation for our chosen triangulation of the doubly periodic domain with n punctures. 
The domain is divided into 2n triangles Uf and Li . Two copies of the domain are shown one above 
the other for clarity. The 3n numbers Xi, yi and Zj are used to count how many crossings a 
lamination makes with each edge of the triangulation. The lamination is assumed to be pulled tight, 
which means that when the curve enters a triangle it must leave through a different edge to the 
one it entered by. Note that this triangulation is not unique, but is chosen for convenience in the 
calculations that follow. 




Fig. 2.3. The lamination inside a triangle T (where T is Ui or Li from Figure \2. can be 
completely reconstructed from the set of crossing numbers {xi,yi, Zi} . To construct the solution, 
each triangle is solved separately, and the triangles are patched together in the obvious way. The 
solution in a single triangle is illustrated. Equation 12.11 determines the number of curves that must 
pass directly between each pair of edges. Up to homotopy, there is only one way to draw these curves 
without crossings. 

infinitely many different ways the domain can be triangulated, so in this paper we 
choose a triangulation that is most convenient for studying the action of braiding 
motions. The notation we use for counting crossings is illustrated in Figure For 
n punctures the resulting encoding {xi,yi,Zi} contains 3n crossing numbers. 

To ensure that each homotopic lamination produces the same set of crossings we insist 
that the lamination is first pulled tight, which means that no loops are allowed where 
the lamination enters and leaves a triangle by the same edge. Under the pulled- 
tight assumption, the set of crossing numbers uniquely identifies the lamination. The 
reason for this is that the path of the lamination can be determined uniquely in 
each triangle, and therefore the global solution is found to be unique by patching the 
triangles together. The unconvinced reader is invited to attempt to find any other 
solution for the given crossing numbers that is not homotopic to the lamination shown 
in Figure E] 

To draw the lamination from the set of crossing numbers we explain the solution 
procedure for each triangle. Consider the triangle T depicted in Figure 1^1 where T 
is Ui or Li from Figure The crossing numbers with each of the three edges are x, 
y and z. In our notation, any part of the lamination passing between edges x and y is 
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counted by T 2 , and, likewise, T y counts lines passing between x and z and T x counts 
lines between y and z. Since we assume the lamination is pulled tight at all times, 
there cannot be any lines that enter and leave by the same edge, so T x + T v + T z 
counts all the lines. A simple bit of bookkeeping shows we must have x — T y + T z , 
y = T x + T z and z = T x + T y , which gives 



rpx 


\{y + z- 


x) 


rpy — 


\ (x + z - 


y) 




\{x + y- 





Given T x , T v and T z , there is then only one way to draw the lines without them 
crossing each other. The T x , T v and T z lines are highlighted using different dash 
patterns for the example in Figure IQ1 The numbers given by (|2.1() are crucial to the 
arguments that follow. 

3. Deformation of laminations under braid operations. In this section we 
describe the effect of braiding operations on a lamination. We start with a lamination 
that is pulled tight, and we record the initial set of crossing numbers {xi,yu Zi}. All 
we have to do now is to determine how these numbers are updated under the action 
of each braid operation. 

To calculate the new set of crossings we simply determine the number of crossings 
with the preimage of each edge. We denote the preimage of edge e by e*, which is a 
curve that becomes e (up to homotopy) after the braid operation. Our argument is 
that the number of crossings of e* before the braid operation has to be equal to the 
number of crossings with e afterwards. This must be the case, as the only way the 
number of crossings could change is if the end of a loop were to cross through the 
edge as it is deformed from e* to e — but by the pulled-tight assumption any such 
loop would have to be wrapped around a puncture, and by construction no punctures 
pass through the edge as it deforms from e* to e. 

It will turn out that we only have to consider explicitly the braid group operations 
Pi, p^ 1 (i = 1, . . . , n) and <7j, a^ 1 (i = 1, . . . , n — 1). We then determine the effect 
of other independent group elements r, and r" 1 by invoking group presentation rules 
and the operation a n (see Figure UTTT b)). This is explained in detail in Section 
for those not familiar with the braid group presentation. For convenience we assume 
in what follows that indices are treated 'modulo' n, so that the puncture to the right 
of puncture n is puncture 1. 

3.1. Crossing update rules for pi. For the braid operation pi, we need only 
consider edges incident on the ith puncture, since the number of crossings with other 
edges will remain unchanged. The relevant preimages are shown in Figure 13.11 as 
dashed lines. Many of the preimages are other edges (e.g. Xi-\ is the preimage y\-\)i 
so the number of crossings is already known. This is not a coincidence: the triangu- 
lation was chosen for this property. The two preimages x*_ 1 and y* require a little 
more attention. In Figure EP1 we illustrate that the number of crossings x*_ 1 can be 
calculated by a simple bookkeeping argument. The number of crossings y* can be 
deduced in a similar way. 
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Fig. 3.1. Preimage illustration for pz- Only edges that are incident with the moving puncture 
can have a change in crossing number. For pz the five edges are X2, V2, xz, yz and 23. In this 
example X2 and 3/3 are the exact preimages 0/1/2 and X3, respectively. The edge 23 is the preimage 
of itself, so in fact its crossing number will not change under pz ■ The preimages of X2 and yz are 
not edges in the triangulation, so more work is required to determine how many crossings are made, 
as shown in Figure T'i. M 




Fig. 3.2. The number of crossings with a preimage that is not part of the triangulation may be 
found by simple bookkeeping. We illustrate here how to determine the number of crossings with the 
X2 preimage. All curves entering X2 and 23 must cross x^ (by the pulled-tight assumption) , unless 
they loop directly from X2 to 23 (shown dotted). The number of these loops is exactly min(f/f , L^). 
Hence the number of preimage crossings is X2 + 23 — 2 miniU^ > -£<! )• Note that curves that cross 
over the preimage and immediately back again are automatically discounted, so that the pulled-tight 
assumption still holds for the updated set of crossing numbers. 



The update rules for all the edges considered are 



x*_ x = x t - X + Zi - 2rnin(f7f_ 1 ,Lf_ 1 ) = x^ x + Zi- X -2min{Uf_ l ,L z l _ l ), 

y*-i =%i-i, 

*i=Vi, (3.1) 
y* = z l + y% -2 wm(U?, Lf) = y t + z l+1 - 2 min([/f , L\), 

* 
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Fig. 3.3. Preimage illustration for 02- The dashed lines show the preimages of all edges whose 
crossing number may change. Edge X2 is its own preimage, but all other preimages require some 
bookkeeping. Figure \3-4\ illu strates how to compute the crossings with x* and j/j* . Computation of 
2/J is done as in Fiaure lH.tA Preimages x^ and y^ are handled in the same way as x^ and y^, by 
invoking rotational symmetry. The difficult preimage problems are and z£, as these preimages 
cross through seven triangles; the solution is described in Section i3.3l 

3.2. Crossing update rules for p^ 1 . We can work out the update rule for 
p~ x by noting the 7r-rotational symmetry of the triangulation about a puncture and 
relabelling the variables in the rules (|3.1|l given above. The update rules for affected 
crossing numbers are 

x *-\ =Vi-l, 

V*-i = Vi-i + Zi - 2min(f7f_ 1 ,£?_ 1 ) - W _i + - 2min(E#_ 1 , LJL X ), 
x* = z t +x t - 2min(Z7f ,Lf) = x t + z l+l - 2 min(Z7? , Lf ), (3.2) 
V* = Xi, 

Zi = Zi. 

3.3. Crossing update rules for <Tj. The same ideas can be used to determine 
the updated crossing numbers following a er^ operation. However, since there are two 
moving punctures, more edges are affected, and the preimages are therefore more 
complicated. An example of a preimage diagram for 02 is shown in Figure 13.31 In 
general, for Oi the affected crossing numbers are aJj-i, Ui-i, Xi, yi, Zi, Xi+i, yi+i 
and Zi+i. This includes the operation <j n which switches the first and last punctures 
across the periodic boundaries (see Figure IrTTT bD : we shall require this operation in 
Sections I3.5H3.6I to find the update rules for and 7^ . 

The edge Xi is its own preimage, and the number of crossings with y* is determined 
using the quadrilateral trick illustrated in Figure 13.21 The number of crossings with 
X i—V Vi—v x i+i an< ^ Vi+i can be found by extending the quadrilateral solution over 
three and four triangles, as illustrated in Figure EOl 

The two remaining crossing numbers to find are z* and z* +1 . Since these preimages 
pass through a total of seven triangles it is non-trivial to determine the number of 
crossings directly. Instead we employ a trick and deduce z* and z* +1 by invoking 
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Fig. 3.4. Two slightly more difficult preimage puzzles for the operation 02- In a similar 
way to in Figure fv. A the number of crossings with preimage xjf is given by the number of curves 
crossing x\ and X2, minus twice the number of loops directly between x\ and X2- The number of 
such loops is exactly min(Lj , US , £f)' Hence x\ = x\ + X2 — 2 min(Lj , U% > L^)- The preimage 
problem for y^ is similar, but involves four triangles and hence a minimum of four numbers, so that 
y*=y 1 + X2-2min(U{,L v 1 ,UZ,L z 2 ). 




Fig. 3.5. The number of crossings with z* and z* +1 are difficult to determine directly from the 
old set of crossing numbers. However they may be calculated using the quadrilateral solution shown 
in Figure Ts. A from the already updated crossing numbers. 



the quadrilateral solution again with the updated crossing numbers. We introduce 
temporary edges p and q directly between punctures i — 1 and i + 1 and between i 
and i + 2, as shown in Figure The preimage of p is yi-i, and the preimage of q 
is Vi+i- Since x*_ x , y*_ lf p, x*, y*, q, x* +1 and y* +1 are known, z* and z* +1 can be 
deduced easily. 

The collected update rules for the affected crossing numbers are 

x*_ x =Xi-i+ Xi - 2 min(LV_ ls Uf , Lf), 
Vt-i = Vi-i +Xi-2 wan(Ul_ 1 ,L v i _ 1 , Uf , L*), 

y*i = Zi+Xi - 2min(t/f ,L 4 2 ) = a* + z l+1 - 2 mm(U t z , Lf), 
z* = <_! + - min(y*_ 1 + - i?, x*_ t + y l - 1 - y*), 

x*+i = x i + - 2 min(E^,if , E7? +1 ), 

= SBi + Vi+i - 2 min(£7f , Lf , C/f +1 , Lf +1 ), 

z *+i = x * + Vi - min(a;* + y i+1 - y* +1 ,y* + y l+ i - x* +1 ). 



3.4. Crossing update rules for crT . Since the triangulation does not have 
a reflection symmetry about a vertical line through the midpoint of two punctures, 
it is not possible to deduce the update rules for a^ 1 by a relabelling in the rules for 
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Fig. 3.6. Preimage illustration for <r 2 . Note that due to asymmetry the preimage problems 
are slightly different to those for 02 ■ 



<7j. An example of a preimage diagram for a^ 1 is shown in Figure l3~Bl The preimage 
curve y* is the most complicated yet, as it passes through ten triangles. However the 
crossing number of all these preimages can be calculated as before using the techniques 
illustrated in Figures IO EHI and 1531 

The update rules for the affected crossing numbers are 



* 

x i-l 


— Xi— 


1 + Xi 


-2mm(U?_ v L*_ 1 ,Uf), 


* 

Vi-l 


= Vi- 


1 + Xi 


-2min(Lf_ 1 ,07)=x i _i+ W -2 


* 


— 






* 

Vi 


= x* 


+ <- 


mm(z* + yi - x*,x* +Vi~ z* +1 ), 


z* 


= Xi 




2mm{Uf,Ly_ l ,U!_ l ,LU,U y i ), 


* 

X i+1 


= Xi 


f x i+ i 


-2min(LlU* +1 ,L* +1 ), 


y'i+i 


= Vi - 


f x l+1 


- 2min(£f , Uf +1 ) = x t + y l+1 - 2 


z i+l 


= Xi 




2mm(LlU? +1 ,Ll +1 ,U? +1 ,L*). 



(3.4) 



3.5. Crossing update rules for ti. In Sections 13.1113.41 we showed how to 
update the set of crossing numbers {xi,yi, z{\ for the braid operations pi, p," 1 , <7j and 
a^ 1 . To complete the set of update rules for any braid we must give the corresponding 
rules for Tj and t" 1 . In performing Tj the ith puncture moves once around the torus 
in the horizontal direction (see Figure fTTTf) . In doing so it passes through many edges 
in our triangulation, so it is difficult to draw the preimages and to derive the number 
of crossings directly. However, we can deduce the update rules for n and r" 1 by 
appealing to group properties that relate <7j, pi and r». In Figure l3~TI we illustrate 
how n is achieved through a sequence of a^ 1 (including a^ 1 defined in Figure lTTlT b)) 
followed by one pj -1 . The other Tj and t~ 1 are produced in a similar manner. A 
computational recipe for each t% operation is given in this section, and t^ 1 in the 
following section. 

To calculate the updated set of crossing numbers {xi, yi, z{\ for Tj do the following: 



10 



M. D. FINN AND J.-L. THIFFEAULT 




Fig. 3.7. The operation n is achieved using a sequence of a. operations, including the 

operation a^ 1 (see Fiaure Tl . il) . followed by a single p^ 1 . This is described in more detail in the 
text. 



1. Use equation H3.4[) to perform, in turn, <J i \ 1 &i-2i • • a i+2 ano - a i+i- Treat 
the indices 'modulo' n, so that erf 1 follows erf 1 . 

2. Relabel xi <— yi «— and 2; <— 2!»+i. This leaves all punctures except 
the ith one in the correct position. 

3. Use equation l|3.2|l to perform pj~ . 

3.6. Crossing update rules for t[~ 1 . To calculate the updated set of crossing 
numbers {xi,yi, z{\ for t~[~ invert the operation of Section 13.51 as follows: 

1. Use equation l|3.1fl to perform pi. 

2. Relabel Xi — > Xi+i, yi — > yj+i and — > 2j+i. This leaves the punctures in 
the wrong position, but the next sequence corrects everything. . . 

3. Use equation 1)3. 3JI to perform, in turn, <7i+i, Oj+2, • • -, Ci— 2 and erj_i. Treat 
the indices 'modulo' n, so that crj follows a n . 

4. Illustrations using simple braids. To illustrate the use of the update rules 
given in Section|31we show how a lamination is deformed under some simple previously 
studied braids. We have implemented the update rules in a short CH — V program, using 
the Gnu Multiple Precision library to allow the number of crossings to grow arbitrarily 
large whilst maintaining exact arithmetic. We also have a Matlab script to draw the 
lamination, one triangle at a time, using the procedure described in Section [3 This 
was used to produce all the figures in this section. In each case our initial lamination 
is a closed loop that passes between the first two punctures, with X\ = y± = 1 and all 
other Xi, yi and Zi set to zero (see the upper- left frame in Figure l4~TJl . 

Figure l4~Tl shows the roll up of the lamination under the repeated action of the planar 
braid o\ with three punctures. In this case it is clear that the third puncture is 
redundant. This braid has zero entropy and is very poor at stirring as it results in 
linear growth of material lines. Note that even though coils form around the pair of 
moving punctures there is always exactly one crossing of the lamination from one copy 
of the domain to the copy above. This is because there was exactly one crossing in 
the initial lamination and this is only a planar braid, so it cannot create any further 
crossings under the pulled-tight assumption. 

Figure FOl shows the result of repeating t-i with three punctures. In this simple braid 
the second puncture moves in a straight line to the left, but catching the lamination 
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Fig. 4.1. Lamination evolution under the braid o~\. This results in linear growth in the number 
of crossings, and zero entropy. 




Fig. 4.2. Lamination evolution under the braid 7"2. This braid also has zero topological entropy, 
but the illustration is provided to verify the correctness of the procedure for performing r using a 
combination of a and p. 

on the puncture exactly once. The illustration is provided to validate the method 
described in Sections I3.5H3.6I for performing r operations using a combination of a 
and p motions. 

The planar pigtail braid c^er^ 1 with three punctures is illustrated in Figure FOl [5]. 
This braid is pseudo-Anosov and has a growth rate per braid letter of |(1 + y/5), 
which is the golden ratio. This 'golden braid' has been proved to have the highest 
topological entropy per braid letter for a planar braid [6]. 

Figure POI illustrates a related cylinder braid oio^o^ 1 17 a 1 with four punctures. This 
is similar to the pigtail braid but wrapped around a cylinder so that the first and 
last punctures are allowed to exchange places. The cylinder braid ctio - ^ 1 with two 
punctures has a growth rate per braid letter given by the silver ratio of 1 + and 
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Fig. 4.3. Lamination evolution under the (planar) golden braid aicr 2 ■ This braid is proved to 
have highest topological entropy per braid letter for a planar braid [6]. 




Fig. 4.4. Lamination evolution under the cylinder braid ai030- 2 cr^ ■ Note that the lamination 
becomes visibly dense more quickly than for the golden braid in Figure [4 . i\ 

this can be proved to be the optimum per braid letter over all cylinder braids. (The 
proof follows that of D'Alessandro [6] using a matrix representation for the cylinder 
braid group.) The silver braid entropy is almost fifty percent higher than the golden 
braid entropy, showing that periodic boundary conditions can be exploited to enhance 
chaos. The extra stretching is clearly visible by comparing Figures f4. 31 and l4~4l 

A final validation of our encoding and update rules is given by checking that any 
lamination is unchanged when it is subjected to the identity braid. For the braid 
group on the torus, one way of writing the identity using a, p and r is a^ 2 p^ 1 t^Pxt^ 1 
which is taken from the group presentation written down by Birman [3]. In Figure l4*31 
the evolution of our test lamination is shown using this braid, and it is unchanged, as 
required. 

5. Calculating topological entropies from laminations. Having completed 
the description of the dynamical system for deformation of a lamination, we now dis- 
cuss how the topological entropy of a braid is related to evolution of a lamination. The 
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Fig. 4.5. Lamination evolution under the toroidal identity braid 2 p 1 1 T2Pit 2 1 ■ The initial 
lamination (not shown) is the vertical line between the first two punctures. The images show, in 

i o 9 1 2 1 2 1 2 1 1 

order, the lamination after a-y , a 1 , a 1 p 1 , a 1 p 1 T2, <r 1 p 1 T2p\ and a 1 p 1 T2pir 2 ■ The 
final lamination is the same as the initial lamination, as required. 

results of Moussafir [16] for the punctured sphere apply here: for an appropriate lam- 
ination, the growth rate of the total number of crossings between the lamination and 
triangulation converges to the topological entropy the braid At. Thus, asymptotically 
we have the estimate 

At = log w + yt + o - log fa + y> + z ^ > f 5 - 1 ) 

i i 

where the starred variables are the updated crossing numbers described in Section [21 
Convergence of X' to the exact braid entropy A appears to be exponential (provided 
the braid has a pseudo-Anosov component), and so in practice very few iterations 
are required to obtain A to double precision. Finite-order braids, which have zero 
topological entropy are detected easily by checking for sub-exponential convergence 
of At. 

As an example, the exact topological entropy of the four-strand braid Uxa^a^ 1 a^ 1 
(see the previous section and Figure 14. 4|) is known to be twice the logarithm of the 
silver ratio A = 2 log(l + a/2) = 1.762747174039086 .... In TableOwe show how the 
total crossing number and the entropy estimate At given by (|5.1() evolve under the 
first few iterations of the braid (using the initial lamination described in Section Q). 
Double precision accuracy of the entropy, which is sufficient for most purposes, is 
reached after about 20 iterations. 

Braids quite naturally arise when considering periodic trajectories of points in a two- 
dimensional flow. For a spatially-periodic flow the space-time plot of these trajectories 
is typically a cylinder or a torus braid, and the entropy of this braid provides a rigorous 
lower bound on the topological entropy of the flow [4,5]. Finn et al. [10] derived such a 
lower bound for the sine flow using Thurston's 'iterate and guess' method to construct 
the train-track for the braid formed by a set of periodic orbits. We will now describe 
how the same result can be found using our lamination approach. 

The sine flow is a time-periodic alternating shear-flow defined on the unit torus 
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iteration 


total crossings 




entropy A^ 




error A — 


1 


24 


2 


,48490664978800 





.72215947574891 


2 


154 


1 


,85889877206568 





.09615159802660 


3 


912 


1 


,77868738766070 





.01594021362162 


4 


5330 


1 


,76546652708556 





.00271935304647 


■5 


31080 


1 


,76321328732169 





.00046611328261 


6 


181162 


1 


,76282713309230 





.00007995905321 


7 


1055904 


1 


,76276089245107 





.00001371841199 


8 


6154274 


1 


,76274952773491 





.00000235369583 


9 


35869752 


1 


,76274757786911 





.00000040383002 


10 


209064250 


1 


,76274724332535 





.00000006928627 


11 


1218515760 


1. 


,76274718592673 





.00000001188765 


12 


7102030322 


1 


,76274717607868 





.00000000203960 


13 


41393666184 


1 


,76274717438903 





.00000000034994 


14 


241259966794 


1 


,76274717409913 





.00000000006004 


15 


1406166134592 


1 


,76274717404939 





.00000000001030 


16 


8195736840770 


1 


,76274717404085 





.00000000000177 


17 


47768254910040 


1 


,76274717403939 





.00000000000030 


18 


278413792619482 


1 


,76274717403914 





.00000000000005 


19 


1622714500806864 


1. 


,76274717403910 





.00000000000001 


20 


9457873212221714 


1 


,76274717403908 





.00000000000000 



Table 5.1 

Convergence of the entropy estimate A^ in 15.11 towards the exact topological entropy A = 
21og(l + y/2) = 1.762747174039086 ... for the silver braid o"i<T3<t 2 cr, . Convergence appears to be 
exponential for a pseudo-Anosov braid, with approximately one extra digit per iteration [16]. 



< x, y < 1. The period is T, with the velocity field given by (sin27ry,0) for < 
t < T/2 and (0,sin27ra:) for T/2 <t<T, with t marking time. This simple flow has 
been well studied because the parameter range < T < 2 gives rich dynamics that 
vary from complete integrability to almost global chaos (with few visible islands in a 
Poincare section) . 

Since the flow is piecewise steady it is easy to construct a map to track the motion of 
points from one period to the next. Consequently, the entropy of the flow for a given 
T can be found quickly and simply by direct numerical simulation of line stretching. 
Once the entropy of the flow is known, it is instructive to see what prediction of the 
entropy is given by considering the braiding of a finite number of particle orbits. 

For general T it is difficult to locate unstable periodic orbits due to the highly chaotic 
nature of the sine mapping; however, for the special parameter T = 1 it is quite easy 
to spot some of the low order orbits. In particular, four period-two orbits are given 
by {(0, ±), (|, I)}, {(i, |), (1, |)}, {(i, 0), (i §)}, and {(f , ±), (f , 1)}. The first pair 
of orbits is depicted in Figure 15. lf al , the second pair in Figure I5.1f b~) . The points 
associated with the first pair of periodic orbits do not move in the second half-period, 
while those associated with the second pair do not move in the first half-period. 

Now we must determine the braid formed by the eight trajectories forming these four 
periodic orbits. We first disambiguate the order of the periodic points by displac- 
ing them slightly along the x axis, as shown in Figure 15.21 Then we encode the 
trajectories in terms of braid group generators, deforming as needed. Deforming is 
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Fig. 5.1. The four periodic orbits considered for the sine flow with T - 
and {(I, |), (1, |)}; (b) {(i,0), (J, |)} and {(§, §), (§, 1)}. 



1. W {(0 A)} 




Fig. 5.2. Illustration of the eight strand braid formed by a set of period two points in the sine 
flow with T = 1. The braid word is <ri cr.7 1 r7/ 1 o~A °~2^ u 7 T 5 u b^ u ^ a 7^ o~2Pqo~q. This 
braid has been shown to have an exact entropy of 1.21875572687 .. . using train-tracks [10]. The 
entropy estimate using the growth of laminations agrees to every decimal place calculated. 



necessary since usually the trajectory does not map directly onto a generator, and 
some intermediate operations must be insterted. For instance, in the second snap- 
shot in Figure IS~21 the generators 03" <7^" er-f 1 are use d to return the point to the 
leftmost position after a t^ 1 operation. After a full period, the resulting braid word 
is (J\02 1 t^ 1 a^ 1 a^ 1 (Ti 1 (JTa^T^a^ 1 ctq 1 cjj 1 p^ 1 a2Pe,(Je,- Using equation l|5.1|l we find 
the entropy of the braid converges to 1.21875572687. . ., which agrees with an alter- 
native calculation using train-tracks [10]. The entropy lower bound given by the braid 
accounts for 82% of the flow entropy of approximately 1.48. 

Since it is difficult to find any exact periodic orbits for general T it is natural to 
ask whether the entropy can be found by considering the braiding of any selection of 
trajectories. Since point motions in the sine flow are piecewise horizontal and vertical, 
and, for the purposes of braiding, can be performed sequentially, it is very easy and 
computationally fast to calculate all the a, p and r operations that occur during each 
half-period of the flow, for an arbitrary number of points. 

For the first half-period, where all motions are horizontal, we record a a operation 
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Fig. 5.3. Convergence of the mean entropy estimate (A') towards the exact entropy A as the 
number of trajectories n in the braid is increased. For large n the error appears to decrease according 
to a power law. For smaller T the convergence is hindered by many of the points being inside periodic 
islands. 

for each change in order of the particle x coordinates. The sign of each crossing is 
determined by the difference in y coordinates for the two points that cross. Special 
attention is required when a point crosses the periodic boundary. Crossing leftwards 
over x = is achieved by t\<j^ x ■ ■ ■ cr^^ this moves the leftmost point one copy of 
the domain to the left, and then, through a sequence of a operations, moves the point 
to position n and undoes all the undesired crossings with the other points. Similarly, 
when the rightmost point crosses over x = 1 the string T^a~_-^ • • - erf 1 is assigned. 
The computation is simpler during the second half-period: all point motions are 
vertical, so no a motions occur. If the ith point in the x direction crosses downwards 
through y = this is labelled pi. Likewise, an upwards crossing of y — 1 is assigned 

P7 1 - 

To estimate the entropy of the flow we evolve a lamination according to the braid 
that results from the point motions. Since the points are not periodic, the growth per 
iteration continues to vary, but the average value of the growth always converges. The 
initial points are chosen at random throughout the domain, and can either live in a 
chaotic region or in a periodic island. Hence the predicted entropy for a given number 
of points may depend on the initial positions. To allow for this we use an estimate 
(A^) obtained by averaging over many realizations with different initial points. 

In Figure l5*31 we show how the braid entropy (A^) converges towards the flow entropy A 
as the number of trajectories in the braid is increased. For the four different values of 
T we have considered, convergence appears to be as a power law for large n. Although 
we are unable to prove that the exact entropy is reached in the limit as n — > oo we 
would expect this since in this limit the braid will contain perfect information about 
the flow. For small values of T (such as T — 0.5 in Figure 15. 3f) the flow contains 
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large islands of regularity which do not contribute greatly to line stretching. In this 
regime n has to become relatively large before there are enough points exploring and 
encoding the dynamics in the small chaotic region. 

Owing to the efficiency with which the braid is determined and analyzed, we are able 
to consider much larger times and number of points n than has been considered in 
previous articles [18]. Another point worth noting is that for T = 1, the mean entropy 
estimate (A^) = 0.45 with n = 8 eight points is much worse than the estimate of 
At = 1.22 found by using the set of eight judiciously chosen periodic points considered 
earlier. This highlights the important role of low order periodic orbits in determining 
much of the nature of the flow. 



6. Discussion. We have derived a dynamical system to compute the evolution 
of a lamination (equivalence class of simple closed curves) under the braiding of an 
arbitrary number of punctures on the torus. The method is essentially a modification 
of the Dynnikov coordinate approach employed by Moussafir [16], but we use a trian- 
gulation encoding that has favorable properties for studying torus braids. Naturally, 
our method also works for the special cases of cylinder and planar braids. However, 
in the planar case our dynamical system still has 3n variables, more than the 2n — 3 
required by Moussafir, so there is clearly some redundancy in this case. Also, the tri- 
angulation of our domain is not unique, so we expect that the details of our method 
are not unique. However our triangulation seems to be the best choice for simplifying 
the arithmetic. 

For completeness we point out that our dynamical system does not work directly 
for n = 2 punctures because in equation (|3.3fl and l|3.4|) this would mean that the 
indices i — 1 and i + 1 refer to the same quantities. There is nothing difficult about 
the case n — 2 though, and in principle one could write down the corrected update 
rules for {xi, y\, z\, x-z, y<z, z-i\. Alternatively, a lazy but convenient workaround is to 
include a redundant third puncture, glued to one of the other two, so that it is slaved 
to its motion. At a little more computational expense this allows one computer code 
to handle all values of n > 2. The case n = 1 is trivial. 

By using similar arguments to Moussafir [16], it can be shown that as the number 
of iterations tends to infinity the logarithm of the number of lamination crossings 
grows at the rate of the braid topological entropy. Though the number of crossings 
grows exponentially fast, our numerical implementation of the dynamical system uses 
a large integer arithmetic library to allow calculation of the entropy to arbitrary 
precision. Convergence appears to be exponentially fast, with approximately one 
digit of accuracy gained per iteration for the braids we have considered. 

In practice, if only a few digits of accuracy are required then double precision floating- 
point arithmetic for the crossing numbers is adequate and can speed up code signifi- 
cantly. A caveat with using floating-point arithmetic is that is destroys reversibility. 
In general by performing a long braid followed by its inverse the initial lamination 
will not be recovered due to exponential growth of small roundoff errors. This is akin 
to irreversibility due to numerical diffusion in trajectory computations in a chaotic 
Stokes flow. However, since the dynamical system itself is exact, the only errors are 
due to roundoff and not to discretization, so even with double precision the dynamical 
system is surprisingly reversible. 
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If finite precision is not acceptable, our method is easily adapted to find exact en- 
tropies. This can be done by shortcircuiting the minimum functions in the update 
rules i|3.1fl - l|3.4fl . Under repeated application of a braid, the dynamical system quickly 
converges to a periodic pattern where it is known in advance which number will be 
taken for each minimum. With this knowledge each call to the min function can be 
replaced by the correct variable and then a linear system can be written down for the 
crossing numbers. The logarithm of the modulus of the largest eigenvalue will give 
the exact entropy. 

This work was motivated by the study of two-dimensional fluid mixing via the braiding 
motion of fluid particle trajectories [11,18]. In this setting we have derived an efficient 
tool that allows practical analysis of large braids. Outside of this particular application 
it is natural to ask whether the method can be generalized for the braid group on 
surfaces of higher genus, as considered by Birman [3]. Since any surface can be 
triangulated, in principle our method could be extended to higher genus. The problem 
anticipated with generalizing the approach is that it is more difficult to visualize the 
preimage problems on a surface with many holes, and also it is not clear in general 
how to exploit group properties to gain a complete set of crossing update rules. 

We thank Colin Cotter and Jacques-Oliver Moussafir for helpful discussions. This 
work was funded by the UK Engineering and Physical Sciences Research Council 
grant GR/S72931/01. 
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